A DENSITY THEOREM ON EVEN FAREY FRACTIONS 



CRISTIAN COBELI AND ALEXANDRU ZAHARESCU 

I Abstract. Let be the Farey sequence of order Q and let and ^^^^^ be the 

, set of those Farey fractions of order Q with odd, respectively even denominators. A 

^ ■ fundamental property of says that the sum of denominators of any pair of neighbor 

C) ' fractions is always greater than Q. This property fails for and for ^^^^ . The local 

density, as Q — > oo, of the normalized pairs {q' /Q, q" /Q), where q' , q" are denominators 
of consecutive fractions in , was computed in |10| . The density increases over 

a series of quadrilateral steps ascending in a harmonic series towards the point (1,1). 
Numerical computations for small values of Q suggest that such a result should rather 
occur in the even case, while in the odd case the distribution of the corresponding points 
, appears to be more uniform. Reconciling with the numerical experiments, in this paper 

we show that, as Q — s- oo, the local densities in the odd and even case coincide. 



> 

(N 

^ ■ 1. Introduction 

m ; 

Questions concerned with Farey sequences have a long history. In some problems, 

in ■ 

O I such as for instance those related to the connection between Farey fractions and Dirichlet 

^ ! constraints. Recently it has been realized that knowledge of the distribution of subsets 
^ ■ of Farey fractions with congruence constraints would also be useful in the study of the 

ix! ■ 

■ periodic two-dimensional Lorentz gas. This is a billiard system on the two-dimensional 
torus with one or more circular regions (scatterers) removed (see j^I], [El, [Z])- Such 
systems were introduced in 1905 by Lorentz [20J to describe the dynamics of electrons 
in metals. A problem raised by Sinai on the distribution of the free path length for this 
billiard system, when small scatterers are placed at integer points and the trajectory of 
the particle starts at the origin, was solved in [H], [0], using techniques developed in ||T], 
[2], to study the local spacing distribution of Farey sequences. 
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The more general case when the trajectory starts at a given point with rational coor- 
dinates is intrinsically connected with the problem of the distribution of Farey fractions 
satisfying congruence constraints. For example, the case when the trajectory starts from 
the center (1/2, 1/2) of the unit square is related to the distribution of Farey fractions 
with odd numerators and denominators. The distribution of the free path length com- 
puted in and ^ is totally different from the one obtained in [7j. This confirms the 
intuition of physicists that, unlike in the case when the trajectory starts at the origin, 
if one averages over the initial position of the particle, the distribution will have a tail. 
It is then reasonable to expect that, in terms of the distribution of Farey fractions, new 
phenomena would be encountered when one replaces the entire sequence of Farey fractions 
by a subsequence defined by congruence constraints. 

As we shall see below, already the case of the subsequence of Farey fractions with even 
denominators presents nontrivial complications. This is mainly due to the fact that in 
there is a large number of tuples of consecutive fractions with odd denominators and 
length growing to infinity with Q. Recently some questions on the distribution of Farey 
fractions with odd denominators have been investigated in |3], ^0], and JH]- It is the 
purpose of this work to derive a result on fractions with even denominators. 

Two fundamental properties of the Farey sequence of order Q state that if a'/ q' < a"/ q" 
are consecutive elements of then a"q' — a'q" = 1, and q' + q" > Q. These properties 
play an essential role in questions concerned with the distribution of Farey fractions. In 
fact, in any problem where one has an element a'/q' of and needs to find the next 
element of call it a" /q", one can use the above two properties in order to determine 
a"/q", as follows. The equality a"q' — a'q" = 1 uniquely determines a" in terms of a',q' 
and q". In order to find q" in terms of a' and q', notice from the above equality that 
a'q" = —1 (mod q'). The inequalities q' + q" > Q and q" < Q show that q" belongs to the 
interval {Q — q', Q], which contains exactly one integer from each residue class modulo q'. 
Only one of these integers satisfies the congruence a'q" = — 1 (mod q'), and this uniquely 
determines q" in terms of q' and a'. Complications arise when one studies a subsequence 
of Farey fractions with denominators in a given residue class modulo an integer number 
d > 2, since in such a case the above two properties fail (see jTUI, [HI for the case of 
fractions with odd denominators). 

In the present paper we study the relative size of consecutive even denominators in 
Farey series. Although the inequality q' + q" > Q fails in this case too, we shall see that 
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Figure 1: From light to dark are represented Figure 2: The sets of type T(3), T(4) and T(r), 
the sets of type T(l). with r > 5. 



the points {q' IQ,q" /Q) have a hmiting distribution inside the unit square [0,1] x [0,1], 
as Q ^ oo. In the following we let 

^?p,even = {^^ 1 < o < ^ < Q, gcd(a,g) = 1, g = (mod 2)} 

and we always assume that the elements of ^q^^^^ are arranged in increasing order. We 
call a Farey fraction odd if its denominator is odd and even if its denominator is even, 
respectively. A new feature that the sequence '^q^^^^ brings in this type of problems 
comes from the following phenomenon. We know from the equality a"q' — a'q" = 1 that 
the denominators q', q" of any two consecutive fractions a' /q' < a" /q" in are relatively 
prime, and in particular not both of them are even. Therefore, if we study the subset 
i^Q odd '-'f '-"^'^ Farey fractions in , we know that any two consecutive elements of 
are either consecutive in , or there is exactly one element of ^^^^ between them. Thus 
there are only two types of situations to consider. By contrast, we may have any number of 
elements from between two consecutive elements of '^^ ^^^^ . This more complicated 

context that arises in the even case, treated in the present paper, forces us to go through 
a significantly larger amount of data than in the odd case. This is also reflected in the 
larger variety of situations that appear in Figures and Tables ^121 below, which show 
various aspects of the distribution in this case. In what follows we study the local density 




Figure 3: The sequences of sets of type T(2). Figure 4: The covering of V{0, 2) by the sets 

of all types. 



of points {q'/Q,q"/Q), with q',q" denominators of consecutive Farey fractions in dq^^cny 
which he around a given point {u, v) in the unit square. We shah show that this local 
density approaches a certain limit g{u,v) as Q — ^ oo, and we provide an explicit formula 
for g{u, v). 

The property of Farey fractions that first drew attention two centuries ago was their 
very uniform distribution of in [0, 1] (see the survey paper 11^ and the references therein). 
It is natural to expect various distribution results, in particular the one obtained in the 
present paper, to continue to hold in subintervals of [0,1]. Let X C [0,1] be a fixed 
subinterval, and denote d'i = dn and = dn H X. We let 

'^Q,even •= {(o'') Q'") ■ q\ denominators of consecutive fractions in d'^^^^^} , 

and consider the set I^Q^even/Q^ which is a subset of the unit square [0, 1] x [0, 1]. If X is 
the complete interval [0, 1], we drop the superscript and write Vq^even = ^^g'lien- 

For each point {u,v) G [0, 1] x [0, 1], we take a small square □ centered at (u, f ) and 
count the number of points from T^'q^evenlQ which lie inside the square □. We shall see 
that, as Q — ^ oo, the proportion of points from T^q ^^^JQ that fall inside □ approaches a 
certain limit. This limit will be proportional to the area of the square □. After dividing 
this limit by Area(n) and letting the side of the square □ tend to 0, we arrive at a limit, 
call it f ), which will only depend on the point (u, f), and possibly on the interval X. 
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Thus, we let 



lim 

q [u,v) := lim ; -— , (1) 

^ ^ ^ Area(n)^o Area(n) ^ ^ 



in which □ C are squares centered at {u,v). We put g{u,v) = g^^'^^{u,v). Theo- 
rem ^ below shows that the limiting local density function g{u,v) exists, and its value is 
calculated explicitly. 

Since the Farey fractions are distributed in [0,1] symmetrically with respect to 1/2, 
the components of the argument of the density g{u, v) will play a symmetric*^ role and we 
shall have g{u,v) = g{v,u), for any {u,v) G [0, 1]^. For convenience, in the following we 
shall use the symbol z for either of the variables m or f and z for the other. Also, to write 
shortly the characteristic function of a system of conditions (equalities or inequalities in 
variables u and v), we denote: 



_ 1, if conditions hold true for {z, z) 

9? (conditions) 

0, else. 



and 



(y9 (conditions) 



1, if conditions hold true for {z,z) = {u,v) and {z^z) = {v,u) 
0, else. 



Theorem 1. The local density in the unit square of points {q'/Q,q"/Q), where q' and 
q" are denominators of neighbor fractions in dq^^cn' approaches a limiting density g as 



"We use tho word symmetric for tuples with components listed in reverse order of one another, and 
also for points situated symmetrically with respect to the first diagonal. 
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Q ^ OO. Moreover, for any real numbers u, v with < u,v < 1, 

1 z + 



g[u,v -- 



if '1 



+H'= = i- (2) 

Figures CMU show how the unit square is covered by countably many polygons, on the 
interior of which the local density function g{u, v) is constant. The explicit value of that 
constant is provided by Theorem^ We remark that for a point {u, v), which is an interior 
point of one of the polygons that form the above covering of the unit square, the second 
sum and the third sum on the right side of (j21) vanish, and in the first sum only finitely 
many terms are nonzero, namely those corresponding to the values of j > 1 for which one 
has simultaneously j < {u + v)/{l — v) and j < {u + v)/{l — u). So, g{u, v) reduces in this 
case to a partial sum of the harmonic series, with more and more terms of the series to be 
counted as the point {u, v) is chosen closer and closer to the upper-right corner (1, 1) of 
the unit square. Additional terms, contained in the second and in the third sum on the 
right side of Q, only appear in the case when {u, v) lies on one of the sides or coincides 
with one of the vertices of one of the polygons that form the covering of the unit square. 

Turning to the same problem on shorter intervals, if X does not have 1/2 as midpoint, 
than the symmetry g-^{u,v) = g^{v,u) is apriori not at all obvious. The next theorem 
provides the stronger result that g-^{u,v) not only exists, for any X C [0, 1], but that it is 
independent of X. 

Theorem 2. Let I C [0, 1] of length > 0. Then 

g^{u, v) = g{v, u), for any {u, v) E [0, 1] x [0, 1]. (3) 

Note that, for any Q > 2, pictures like those from Figures A and B are symmetric, 
point by point, with respect to the first diagonal, and this also holds in the case of 
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and d'^^^^^, for any interval X symmetric with respect to 1/2. If X is not centered at 
1/2, then the corresponding pictures are never symmetric point by point, but Theorem|21 
guarantees that their hmits, as Q ^ oo, are symmetric. 

By comparing Theorem Q above with Theorem Q of flO\, we find that the hmiting local 
density function g{u,v) is the same for both subsequences dg^dd dq^^^n- Therefore, 
Corollary 1 and Proposition 1 from ^U] also hold for ^^^^ . 

Let X'(0,2) denote the set of limit points of sequences of pairs i{<ln/Qn,(ln/Qn))n£Ny 
with Q„ — > oo and g^, g" denominators of consecutive fractions in dg^^^n- 

Corollary 1. The set V{0, 2) coincides with the quadrilateral bounded by the lines: y = 1, 
X = 1, 2x + y = 1, and 2y + x = 1. 

We remark that the more uniform distribution of X'socodd compared to X'soo.even in 
Figures A and B is due to a combination of two factors that occur for Q small. The first 
one is the preponderance in of pairs of type T(0), that is, of neighbor odd fractions 

in , and the second is the fact that the cardinality of is about twice as large as 

the cardinality of dg^^^n- 

The next corollary provides the probability that a neighbor pair of denominators in 
^„ is small, in the sense that their sum does not exceed Q. 

Corollary 2. The probability that the sum of neighbor denominators of fractions from 
i^Qcvcn is < Q approaches 1/6, as Q ^ oo. 

In order to prove Theorem^ we begin by presenting in Sections El and El some geometric 
prerequisites, in particular the tessellation of the Farey triangle, and then continue in 
Sections |^ and El with the study of different types of pairs of neighbor denominators of 
fractions in dg^^^n- These conclude with Theorem El in which gr{u,v), the local density 
at level r, is written as a sum of a series of special quantities assigned to the pieces of the 
tessellations of the Farey triangle. This result plays a key role in the proof of Theorem Q 
from Section ini where we put together all the pieces. Remarkably, these pieces fit into a 
large mosaic composed by a series of superimposed puzzles at odd levels (see Figures EHHI) 
recovering in the even case the same density as that obtained in JUj in the odd case. 
In the last section we employ a technique which makes use of estimates for Kloosterman 
sums [ini , [221 , in order to prove Theorem |21 




2. Some Geometry of the Farey Fractions 

We state here the basic properties of the Farey series needed for the rest of the paper. 
The first one, already mentioned in the Introduction, says that if a' /q' < a" /q" are 
neighbor fractions in 3"^, then 

a"q' - a'q" = 1 . (4) 

Next, suppose that a' /q' < a" /q" < a'" jq"' are consecutive elements of 5^^. Then, the 
middle fraction, called the mediant, is given by 



a" a' + a'" 



(5) 



q" q' + q'" 

This shows that the mediant fraction is reduced by an integer k, called the index of the 
Farey fraction a'/g', that satisfies: 

'Q + q' 



k = — = — = a q -aq 



q- 

We state a third important property in the following lemma 



q" 



(6) 



Lemma 1. The positive integers q',q" are denominators of neighbor fractions in if 
and only if {q',q") G Tq and gcd{q' , q") = 1. Also, the pair {q',q") appears exactly once 
as a pair of denominators of consecutive Farey fractions. 
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For the proof of relations (jH) and®, observed for the first time in particular cases by 
Haros and Farey, we refer to Hardy and Wright JH] 5 while for Lemma ^ and further 
developments, see Hall [H], [THj, Hall and Tenenbaum jT7j, and Hall and Shiu [TB] . 

Now let (g', g", g'", . . . , q^^^) denote a generic /i-tuple of denominators of neighbor frac- 
tions in Then, we see that relation Q can be employed together with Lemma Q 
to obtain a characterization of any such /i-tuple in terms of its first two components. 
Moreover, while any pair (g', g") with coprime components < Q does appear exactly once 
as a pair of neighbor denominators of Farey fractions, the components of longer tuples 
must satisfy additional conditions in order to appear together as neighbor denominators 
of fractions in 3"^. We write these conditions using the index. 

For any positive integer k, we consider the convex polygon defined by 

r^, := {{x,y): 0<x,y<Q, x + y>Q, ky < Q + x < {k + l)y} . 

These are quadrilaterals, except for k = 1, when 7^ ^ is a triangle. The vertices of 
7^^ are (0,(5), ^) ; (Q^Q), and for any k > 2, the vertices oiT^^ are (Q, ^) ; 

; B2) 5 S)- Scahng by a factor of Q, for any > 1 we get 

Tk = Tq^/Q, bounded polygon inside the unit square, which is independent of Q. 

On each the index function, defined by (x, y) i— > [^^] , is locally constant. Therefore 
the polygons 7^ can be used to describe the triplets of neighbor denominators of Farey 
fractions. We state this lemma. 

Lemma 2. The positive integers g', g", q'" are denominators of consecutive fractions in 
and k = if and only if (g', q") € 7^ j. and gcd(g', q") = 1. 

We remark that the sets 7^, with k > 1, are disjoint and they form a partition of T. 

In many instances, one can estimate the number of Farey fractions with a certain 
property by counting the number of lattice points in a suitable domain. In this respect, 
the following lemma, which is a variation of Lemma 2 from [1], is very useful. For any 
domain C we denote: 

A^odd,odd(^^) ■■= {{x,y) e finZ"^ : X odd, y odd, gcd{x,y) = l} , 

A^even,odd(fi) := # {ix,y) E Vt Hi? : X even, y odd, gcd(x,y) = l} , 

A'odd,even(f^) := {{x , y) E Vt f] : X odd, y even, gcd(x, y) = l] . 
These numbers can be estimated by using Mobius summation. One has the following 
asymptotic formulas. 
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Lemma 3 C^, Corollary 3.2). Let Ri, R2 > 0, and R > mm{Ri, R2) ■ Then, for any 
region Q C [0, Ri] x [0, R2] with rectifiable boundary, we have: 

NoddMdi^) = 2Area(fi)/7r2 + 0(CR,f,), 
iVodd,even(f^) = 2Area{n)/n^ + 0iCR,n), 
where Cr^q = Area(n)/i? + R + length((91]) hgR. 

Then, for Q as in Lemma El one immediately gets 

iVeven,odd(^^) = 2Area(^]) /tt^ + OiCn,n)- 
3. The polygons T 

k 

We consider the map T : T ^ T, defined by 



l+x 



y-x , 



y 

which was introduced and studied in j3] . This transformation is invertible and its inverse 
is given by 



T-\x,y) 



l+y 



X — y,x 



One should notice that if a'/q' < a" /q" < a'" /q'" are consecutive elements in "^q, then 
T{q'/Q, q"/Q) = iq"/Q, q"'/Q). Then, for any k = {h, . . . .K) E we put 

Tk = n n ■ ■ ■ n T-'^'%^. 

We use the notational convention of dropping extra parentheses, so for example when 
r = 1 and k G N*, we have T(^k) = ^ = {{.x,y) G T: \}-^'\ = Also, to avoid double 
subscripts, at small levels r we often write k,l, . . . , rather than ki, k2, ■ ■ ■ 

We remark that at any level r, the convex polygons 7k are pieces of a partition of 
T. The structure of these partitions is studied in jT^j. Here we only need the polygons 
assigned to tuples k, whose components have special parities. These are those tuples k 
with all components even, except for the first and the last when r > 2. We call these 
tuples admissible, and for each level r, we denote by A{r) the set of admissible r-tuples. 

The polygons needed in the sequel are listed in Table Q Notice that here 7^ is a 
quadrilateral, except for k = (1, 2, 4, 1) and k = (1, 4, 2, 1), when it is a triangle. 
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Table 1: The polygons 7k, for all admissible k. 



Level 


k 


No. of vertices 


Vertices of 7k 


1 


fc > 2, even 


4 




2 


1,3 


4 


(i,|); (|,f); (1,1); (1,1) 


2 


(1,0, with I > 5 odd 


4 


/i — 3 '-IV (1 — 2 I ■\. (l-l i\. /I — 2 -1-. 
U + l'i+l/' U+2'!+2^' Vi+l'-"-/' i ' 


2 


3,1 


4 


(1,1); (l,f); (1,1); (1,2) 


2 


{k, 1), with fe > 5 odd 


4 


fi£^ (_k 2_Y Q Q 2-) 


3 


1,2,3 


4 


(i,f); (1,1); (1,1); (1,1) 


3 


3,2,1 


4 


(1,1); (1,2); (1,4); (1,1) 


3 


1,4,1 


4 


(757)5 (3' 3)' (75-^)' (2'-^) 


3 


(1, 1, 1), with I > 6, even 


4 


(1-3 (1-2 I \. (l-l -ix (1-2 i\ 


4 


1,2,2,3 


4 


(959)5 (757)5 (551)5 (75 1) 


4 


3,2,2,1 


4 


5 ' 5/' V7' T/' V ' 9/' ^ ' 


4 


1,2,4,1 


3 


(it); (ii); (f,i) 


4 


1,4,2,1 


3 


(1,2); (1,1); (1,1) 


r 


1,2,. ..,2,3 


4 


2r+l ' 2r+l /' V 2i — 1 ' 2i — l/' V 2i — 3' /' V2t--1' / 


r 


3,2,. ..,2,1 


4 


/2r— 5 r— 2 v f2r—3 r— 1 \. /"-i r+1 y /l ^ 
V2T--3' 2T--3-'' V2t--1'2i — 1/' V""-, 2t-+1 z' V-'-i2t--1'' 



We shall denote Tq^^ :— Q%^. 

4. The components of X>(0, 2) 

Since any two consecutive denominators of fractions in are coprime, it follows that 
between any two neighbor fractions in there must be at least one odd fraction from 

We note that the number of these odd fractions may be quite large as Q increases. We 
classify the pairs of neighbor even fractions according to the number of odd intermediate 
fractions. So, we say that the pair (7', 7") of consecutive fractions in ^^^^ is of type 
T(r) if there exist exactly r odd fractions between 7' and 7" in ^q. In this case, we 
shall also say that the pairs (g', q") and {q'/Q, q"/Q) are of type T(r), where q', q" are the 
denominators of 7', 7", respectively. 

Geometry inside the Farey triangle helps one locate the points corresponding to such 
pairs. Next, we determine, one by one, the contribution of pairs of each type T(r) to 
T>{0, 2). We denote by £{r) the contribution to T>{0, 2) of points of order T(r). Then, we 
have 

V{0,2) = \JS{r). (7) 

r=l 

In the following we find exphcitly each set £{r). 
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4.1. Points of type T(l). By it follows that pairs of fractions of this type have 
as denominators the end points of a triple {q\q",kq" — q'), with q' even, q" odd and 
ki = k = [^77^] even. This means that, for any even k > 2, we need to retain the lattice 
points in the domain 

UQ,k={{q',kq"-q'): (g', g") ^ T^, J , 

with q' even, q" odd, and gcd(g', q") = 1. In the limit, when Q oo, these points produce 
a subset of V^O, 2) that is dense in the quadrilateral with vertices 

Thus, we have 

oo 

S{1)= U U,.^ (8) 

k=2 
k even 

We remark that all Uk, not only those with k even, contribute to T>{1, 2) (see [TU]). 

4.2. Points of type T(2). These points come from 4-tuples of parity (e, o, o, e) and this 
requires both k and / to be odd. Let 

Wq,., = { (g', l{kq" - q') - q") : (g', q") E T^,,, } , 

and we need to pick up the lattice points in UQ^k,i where g' even, q" odd, and gcd(g', g") = 1. 
In the limit, when Q oo, these points produce sequences of subsets of V^O, 2) that are 
dense in each of the quadrilaterals with vertices 

Wi,3 ={(1/5,1); (2/7,4/7); (1/2,1/2); (1/3,1)}, 

W3,i={(l/2,l/2); (4/7,2/7); (1,1/5); (1,1/3)}, 

^-=U^>i); f^>A); Ct^^'t^I f^>l)l for/>5, 

• IV/ + 1' /' V/ + 2'/ + 2/' V/ + l'/ + l/' V Z ' //' - ' 



k + V k + lJ' \k + 2' k + 2J' \ ' k+lJ' \ ' k 
Then 

oo oo 

S{2) = Wi,3 U ^3,1 U U U U Ui,i . (9) 

A:=5 «=5 
k odd / odd 



'^Most often, we denote a polygon by its vertices and, depending on the context, the same notation 
will be used to indicate its topological closure. 
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4.3. Points of type T(3). These points are produced by 5-tupIes of parity (e, o, o, o, e) 
and this requires both k and m to be odd and / even. Let 



Uq,,,,^ = {{q',m{l{kq" - q') - q") - (kq" - q')) : (g', g") G T^^^^.J , 

and we need to pick up the lattice points in UQ^k,i,m where q' even, q" odd, and gcd(g', q") = 
1. In the hmit, when Q ^ oo, these points produce the sequence of subsets of 1^(0,2) 
that are dense in each of the quadrilaterals with vertices 

ZYi,2,3={(l/r,l); (1/5,3/5); (2/7,4/7); (1/5,1)}, 
W3,2,i ={(4/7,2/7); (3/5,1/5); (1,1/7); (1,1/5)}, 
ZYi,4,i ={(2/7,4/7); (1/3,1/3); (4/7,2/7); (1/2,1/2)}, 

"^'''^={(7TI'7Tl)' (ttI'ItI)' (iTI'm)' (^'^)}' 

Then 

oo 

£(3) =Wi,2,3UW3,2,iU U (10) 

I even 

4.4. Points of type T(4). These points are produced by 6-tuples of parity (e, o, o, o, o, e) 
and this requires both k, n to be odd and m, I even. Let 

l^Q,k,i,m,n = { {q, n[m[l{kq" - q) - q") - {kq" - q)) - [l{kq" - q) - q") : 

(9',/)e V,_}. 

As before, we need to pick up the lattice points in lAQ^k,i,m,n with q' even, q" odd, and 
gcd(g', g") = 1. By Tabled we know that there are only four such sets. In the limit, as 
Q — > oo, we obtain two quadrilaterals and two triangles: 

Wi,2,2,3 ={(1/9,1); (1/7,5/7); (1/5,3/5); (1/7,1)}, 
^3,2,2,1 ={(3/5, 1/5); (5/7,1/7); (1,1/9); (1,1/7)}, 
Wi,2,4,i ={(1/5,3/5); (1/3,1/3); (2/7,4/7)}, 
Wi,4,2,i={(l/3,l/3); (3/5,1/5); (4/7,2/7)}. 

Then 

^(4) = Wi,2,2,3 U ^3,2,2,1 U Wi,2,4,l U Wl,4,2,l ■ (H) 
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4.5. Points of type T(r), r > 5. These points are produced by the (r + 2)-tuples 
(g', . . . , g^*""^^)) of parity (e, o, . . . , o, e). It turns out that the only possible corresponding 
r-tuples k are (1, 2, ... , 2, 3) and its symmetric (3, 2, ... , 2, 1), with (r — 2) 2's between 
the endpoints. This allows us to find a closed formula for q^'^'^'^\ We find that = 
— (2r — \)q' + 2q" in the first case, and = — g' + 2g" in the second case. Thus, we 

define 

Wq,i,2,...,2,3= {(g,-(2r-l)g' + 2g"): ^ ^,,...,,,3} , 

and 

^^Q,3,2,...,2,i = { (g', -g' + 2g") : (g', g") G r^,3..,.....,, } • 

As before, only the lattice points with q' even, q" odd, and gcd(g', g") = 1 should be 
considered, and in the limit, when Q 00, for any given r, we obtain the quadrilaterals: 



1 \ / 1 2r-3\ / 1 2r-5 



2r-5 1 \ /2r-3 1 \ 1 \ 1 



Then 

^(r) = Wi,2,...,2,3 U W3,2,...,2,i , for r > 5. (12) 

In Figure ini one can see a representation of T>{0,2) covered by U^:. In addition, the 
union from the right hand side of ((Tj) gives a first hint on the local densities on V{0,2). 
Complete calculations are postponed to Section El 

5. The density of points of type T(r) 

Let g{x, y) be the function that gives the local density of points {q' /Q, 1" IQ) in the unit 
square as Q — > cxo, and let gr{x,y) be the local density in the unit square of the points 
{q'/Q,q"/Q) of type T(r), as Q ^ 00. At any point {u,v) G (0, 1)^, this local density is 
defined by 

#(DnpQ(0,2)/Q) 
o"So #»q(0,2) 

qJu.v) := lim -— , (13) 

' ^ Area{n)-.o Area(n) ^ ' 

where □ C are squares centered at (u, v). They are related by 

00 

g{u,v) = ^gr{u,v) , (14) 

r=l 
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provided we show that each local density gr{u,v) exists as Q ^ oo. In the following we 
find each gr{u, v). The proof generalizes that of Theorem 1 from piU Section 3.2]. 

5.1. Generalities on Qj.. Let r > 1, let {xo,yo) be a fixed point inside the unit square 
(0, 1)^, and fix a small rj > 0. In the following, the superscript C indicates the last 
element of a tuple. For example, = Qril'il") is the denominator of the (r + 2)- 
nd fraction in 5^^, starting with g', g". We denote the square centered at (xo,?/o) by 
□ = □^(xo,2/o) = {xq -ri,XQ + vi)x {yo -r],yo + r]). 

By definition, any pair (g', q^) of type T(r) is generated by an (r+2)-tuple (g', q", . . . , g^) 
of denominators of consecutive fractions in with g', g^ even, and the rest of the 
components odd. We consider the set Bq of pairs (g', q") of type T(r) for which the 
corresponding point {q' /Q,q'^IQ) falls in □, that is, 

g' even, q" odd; k(g',g") G ^(r), {q,q^{r)) G Q ■ □ 

The cardinality of Bq is #;B<3 = iVeven,odd(^Q), where Qq = ^^(r) = QQ{xo,yo,r]){r) is 
given by 

k(x, ?/) G ^(r), (x, (x, y)) G Q ■ □ 

Here x^ = x^{x,y) is the (r + 2)-nd element of the sequence defined recursively by: 
x_i = X, Xq = y and Xj = kjXj-i — Xj-2, where kj = [-^r^] , for 1 < j < r and 
k{x,y) = {ki, ...,kr)e A{r). 

Multiplying by l/Q, we obtain the bounded set 

n{r) = {{x,y) G (0,1)2; 3. + ^ > 1^ ^^^y^ ^ ^(^)^ G □} , 

and Q ■ Q{r) = ^^(r). We are interested in the area of fi(r), since, by LemmaEl we know 
that 

, , 2Q2Area(fi(r)) , , 
#Sq r = ^ + 0{Q logQ . 15 

Next we split Q{r) into the pieces given by the shadows left on it by each of the sets 



Wk(r) := G (0,1)2; k^(^^2/) = k}^ for k G ^(r) 
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Since n{r) C T and Ui,{r) n T = Tk, it follows that fi(r) n U^{r) = %n P^{r), where 
^^{r) = Vk^r{xo, 1/0,1]) is the parallelogram 

Vk{r) = {{x,y) e M^: {x,x^{x,y)) G □^(xo,l/o)} • 
Thus, we have obtained 

Area(fi(r)) = ^ Area(7; n Pk(r)) . (16) 

A compactness argument shows that, although A{r) may be infinite, only finitely many 
terms of the series are non-zero. Our next objective is to make explicit their size in terms 
of the position of (xo,?/o)- But first we need to establish a concrete expression for Vk{r). 

5.2. The index Pr.(k) cind the parallelogram Vk{r). First, we define the sequence of 
polynomials Pr(k) by: Po{-) = 1, Pi{ki) = fci, and then recursively, for any r > 2, 

Prih, . ..,kr) = krPr-l{ki, K-l) - Pr-2{ki, . . . , kr-2)- (17) 

The first polynomials are: 

P2(k) = kik2 - 1; 

Ps{k) = kik2ks -ki- ks; 

P4(k) = ^1/^2^3/^4 — kik2 — kiki — k^ki + 1; 

P5(k) = kik2ksk4kr, - ^1/^2^3 - kik2kr, - kik^k^ - k^^k^k^ + fci + /C3 + k^. 

Some other particular values we need later are: 

p,(l,2,...,2) =1, forr>2, 

p,(l,2,...,2,3) = 2, forr>2, (18) 

Pr(2, ...,2,3) = 2r + l, for r > 2, 
Also, we remark the symmetry property: 

Pr{kr, ...,ki)= Pr{ki, . . . , kr) . (19) 

More on this fundamental sequence of polynomials can be found in |12j . 
Next, let us notice that x^{x,y) is a linear combination of x and y: 

x^{x, y) = Pr{ki, . . . , kr)y - Pr-i{k2, . . . , kr)x . (20) 
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Returning now to our parallelogram, by (j^UI) we see that this is the set of points 
{x, y) G that satisfy the conditions: 

Xq — Tj < X < Xo + Tj, 

yo-r] < pr{ki, kr)y - Pr-i{k2, kr)x <yo + r]. 
From this, we see that the area of V]i{r) is 

Area(a(r)) = ^ , (21) 

and its center has coordinates 

^ / N f Pr-i{k2,...,kr) 1 \ 

^'^^''^ = 1'^°' 'VJk. — kV'"' + V(k. — J • ^^^^ 

5.3. The density gr{u,v) II. In order to obtain a concrete expression for the density, 
one requires an explicit form of the series in (fT^ . Since rj > can be chosen as small 
as we please, it follows that the summands there depend on the position of {xo,yo) with 
respect to 7^. In the following, we shall assume that rj is small enough. We may also 
assume that k is bounded, since all the parallelograms Vk{r) are contained in the vertical 
strip given by the inequalities xq — t] < x < xq + f], and since only finitely many polygons 
7k intersect this strip. 

Let now k be an admissible r-tuple. We check what happens when C]f_{r) lies on the 

o 

interior 7k of 7k, on the edges of 7k, or in the set Vik) of vertices of 7k. 

Firstly, if Ck{r) e Tk, it follows that Vk{r) C Tk, so Area(Tk n Vk{r)) = Area(Pk(r)). 
Secondly, if Ck{r) G 97k \ V(Tk), then Area(7k fl Vk{r)) = Area('Pk('"))/2, since any line 
that crosses V]i{r) through its center cuts the parallelogram into two pieces of equal area. 
Thirdly, suppose Ck{r) G V^(7k). Then Area(7k fl Vk{r)) depends on the angle formed 
by the corresponding edges, and these angles may differ for different vertices of 7k or for 
different values of k. We have collected all the results in Table |21 In the calculations, we 
have made use of the relations (fTHj) . (fT^ and (plj) . For each V G ^(7k) and Ck(r) = V, 
for a concise presentation, we have translated 7k fl Vi(_{r) with a vector V to the origin. 
Also, we remark that the size of 7k fl V]i{r) is always proportional to 77. It follows that 
C(k(y) '■= Area(7k fl V\^{r))/rf is independent of rj. 

The nice thing about this calculation is that it has an error-correcting check. This is 
due to another remarkable property of the parallelogram, which implies that the entries 
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in the column of ak{V) satisfy the relations 

y a^iV) = I ^^^^^^ = is a quadrilateral; ^^^^ 

\ ^^^^ = if is a triangle, 

for each admissible k. Here the summation is over all the vertices of 7k. For example, if 
k = (1, /, 1), with / > 6 even, we have 

2/-1 1 + 2 2/-1 I 4 

+ T-. ^ + -7^. TTT + 



2{l-l)l {1-2)1 2(1-1)1 {l-2)(l-l) 1-2 
and for k = (1, 2, 4, 1), we have 

7 3 19 _ ^ 
24 ^ 40 ^ 30 ~ 



Table 2: The vertices and the area of the polygons 7k n -Pk(?^') for Ck(r) equal to each of the vertices 
of 7k. Each polygon was scaled by l/?7 and translated into the origin. 



k 


\'(-rt(-x \ ' = f 'k 


Vortices of i (Tk n i'k( /■)) - r'k( r) 


<>k(\-) 


k>2, even 




(0,0); (^,-^); (1,0); (l,i) 


2fc(fc+l) 




(Jl 2--) 

\k+2' k+2l 


(0,0); (1, iqrj); (1, fc); (— ■fcq:T,iq:T) 


fc+2 
fc(fc+l) 




(i>4t) 


(0,0); (0,1); (-1,0); (-1,-^) 


2fc+l 
2fc(fc+l) 




(I'D 


(0,0); (-1,-i); (-1,-1); (0,-i) 


1 
fc 


1,3 


(1 ±) 


(0,0); (i,-i); (1,1); (1, f ) 


9 

20 




(1 5^ 

W' 7' 


(0,0); (1,|); (1,2); (-1,1) 


19 
30 






(0,0); (-1,0); (-1,-1); (-1,-|) 


1 

3 






(0,0); (-1,-1); (-1,-2); (1,0) 


7 
12 


{1,1), with / > 5 odd 




(0,0); (7^,7^); (1,1); (1,1^) 


I 

(i-l)(i+l) 




(1-2 I \ 

y 1+2 ' 1+2 ' 


(0,0); (l,i±l); (l,i±i); (-j^,-,^) 


2i^+5! + l 
2((-l)(((+l) 






(0,0); (-1,0); (-1,-1); (-1,-^) 


1 
I 






(0,0); (-1,-T^); (-1,-l^i); (i,0) 


21+1 

21(1-1) 


3,1 


^2' 2) 


(0,0); (1,- 


t); (1,0); (1,1) 


1 
3 






(0,0); (1,2) 


(1,1); i-hl) 


19 

30 




(1,1) 


(0,0); (0,1) 


(-1,0); (-1,-1) 


9 

20 




(14) 


(0,0); (-1,-1); (-1,-1); (0,-i) 


7 
12 


(fc, 1), with A: > 5 odd 


/fc-l 2 \ 
V + 1 ' k + 1 1 


(0,0); (1,-1); (1,0); (1,1) 


1 

k 




(— —) 

V fc+2 ' fe+2-' 


(0,0); (1,^); (1,-^); (-1,1) 


2fe^+5fe+l 
2(fc.-l)fc(fe+l) 




(1'4t) 


(0,0); (0,-^); (-1,0); (-1,-^) 


k 

(fc_l)(fc+l) 




(1-1) 


(0,0); (-1,-1); (-1,-^); (0,-^) 


2fc+l 
2fc(fc+l) 


1,2,3 




(0,0); (1,-1); (1,2); (1, |) 


13 

28 




(hi) 


(OA)): (Ll); (i.:!): (-4.4) 


1:! 
21 



continued on next page 
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k 


Vertex V = Ck(r) 


Vertices of i(Tk n Pu{r)) - a(r) 


ak(V^) 




(1,1) 


(0,0); (-i,0); (-1,-2); (-1,-|) 


11 

30 




(il) 


(0,0); (-1,-1); (-1,-3); (i,0) 


11 
20 


3,2,1 


(- 2) 


(0,0); (|,-i); (1,0); (1,2) 


11 
30 




(3 2-) 
V 6 ' 6 / 


(0,0); (1,1); (1,1); (-i, i) 


13 
21 




(1,1) 


(0,0); (0,i); (-1,0); (-1,-f) 


13 
28 




(1,1) 


(0,0); (-1,-1); (-1,-1); (0,-i) 


11 
20 


1,4,1 


(14) 


(0,0); (i,-i); (1,1); (1, |) 


11 
30 




(if) 


(0,0); (1,2); (1,2); (-1,1) 


3 




(1,1) 


(0,0); (-1,0); (-1,-1); (-1,-1) 


11 
30 




(il) 


(0,0); (-1,-1); (-1,-2); (1,0) 


2 
3 


(1,/, 1), with I > 6 even 




(0,0); (1,1): (1,1): (1,^) 


2i-l 

2(i — 1)! 




V Z4-2 ' Z+2 / 


(0,0); (1,^); (1,^); (-f,-i) 


1+2 
(1 — 2)1 




(l+f,!) 


(0,0); (-,^,0); (-1,-1); (-1,-^) 


21-1 
2(1-1)1 




(¥,1) 


(0,0); (-1,-T^); (-i,-i42); (7^,0) 


I 

(l-2)(l-l) 


1,2,2,3 


fl ») 


(0,0); (1,-1); (1,3); (1,1) 


17 
36 




W, 7' 


(0,0); (i,|); (1,3); (-1,1) 


65 
72 




(il) 


(0,0); (-1,0); (-1,-1) 


5 

56 




(1,1) 


(0,0); (-1,-1); (-1,-4); (1,0) 


15 
28 


1,2,4,1 


(1 4-) 
V6, 6/ 


(0,0); (1,1); (1,2); (1, |) 


7 

24 




(1,1) 


(0,0); (-1,0); (-1,-1) 


3 

40 




(f,l) 


(0,0); (-1,-1); (-1,-3); (1,0) 


19 

30 


i, 1,2, i 


(it) 


(0^0): (Ll): {Ll) 


40 




(1,1) 


(0,0); (-1,0); (-1,-1); (-1,-1) 


7 

24 




(1,1) 


(0,0); (-1,-1); (-1,-2); (1,0) 


19 
30 


3,2,2,1 


(3 

V 5 ' 5 / 


(0,0); (1,1); (l,f) 


5 

56 




(S 3N 


(0,0); (1,|); (1,1); (-1,0); (-1,-1) 


65 
72 




(1,1) 


(0,0); (0,1); (-1,0); (-1,-|) 


17 
36 




(1,1) 


(0,0); (-1.-I); (-1,-1); (0,-i) 


15 
28 


1,2,. ..,2,3 


2l_) 

V 2r+l ' 2r+l / 


(0,0); (^,-^); (l,r-l): (l/^) 


4r + l 
4(2r + l) 




/ 1 2(r-l)x 
V2r— 1' 2r — 1 / 


(0,0); (1,2^^); (l,r-l); (l,r); {-^^,^^) 


14r+9 
8(2r + l) 




(27^,1) 


(0,0); (-27^,0); (-1,-2:^) 


2r-3 
8(2r-l) 




(27^,1) 


(0,0); (-1,-21^); (-1,-r); (^,0) 


4r-l 
4{2r-l) 


3,2,. ..,2,1 


/ 2r — 5 r — 2 \ 
V2r-3' 2r-;V 


(0,0); (1,1); (1,^) 


2r-3 
8(2r-l) 




/ 2r — 3 r — 1 \ 
V 2r-l ' 2r-l 


(0,0); (1,2;^); (1,1); (-1,0); (-1,-1) 


14T-+9 
8C2T-+1) 




fl '• + 1 ) 


(0,0); (0,1); (-1,0); (-l,-^;^) 


4r-+l 

4(2r + l) 




(1-2;^) 


(0,0); (-1,-^); (-1,-1); (0,-1) 


4r-l 
4(2, — 1) 
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These observations combined with ()21|) put (fT^ into the form 

Area(fi(r)) = V E ^ + E ^ + E ""^W ' (24) 

We now immediately obtain a corresponding expression for An apphcation of 
Lemma El (see also Lemma 0] below) provides 

#i?g,_=QVvr' + 0(QlogQ). 

Using (|15|) and the fact that the number of points {q' /Q-, q" /Q) from (0, 1)^, where q' and 
q" are denominators of two consecutive elements from 5„ is — 1, we have: 

^ ^Q,cvcn ^QjGvcn ' 

/ / Qrix, y) dxdy = hm ,f ^^^""^ ^ = 2Area(n(r)) . (25) 

Then, by the Lesbegue differentiation theorem, we have gr{xo, yo) = lim^^o Area(i7(r))/ (4?7^) 
which combined with gives the following result. 



Theorem 3. For {xo,yo) G [0, 1]^ and any integer r > 1, we have: 

9r{xo,yo)= E ^ + 2E "k(r), (26) 

where the sums run over r— tuples k which are admissible. 

We call the generic term in the first sum the kernel of k. Notice that if C]f^{r) is on the 
boundary the terms added in the second sum are equal to half of the kernel, and in the 
third sum the kernel distributes in sizes proportional to the angles of different vertices of 
the polygon 

6. Proof of Theorem [T] 

Our starting point in the proof of Theorem ^ is formula ()26p . in which we need to 
make explicit the conditions of summation in terms of the variables. Since some of the 
polygons with all components of k small do not follow the general pattern, at this point 
we can not get an explicit closed formula for the density, even if we restrict to pairs of 
a given order. Consequently, we further split the sums on the right hand side of (j2EI)- 
For this, we collect in h^{xQ,yo) the sum of the terms from the right hand side of 
with r = 1 and k = > 4, even. Similarly, for r = 2, h^{xQ.,yQ) will be the sum over 
k = (1, Z) and k = (/c, 1) with k,l > 5, both odd, while for r = 3, h^{xo,yo) will be the 
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sum over k = (1, /, 1), with / > 6, even. The remaining terms will be collected separately 
in hf{xo,yo), h2{xQ,yQ) and /if (xq, 2/o)- Thus, we have: 



gj{xo, Vo) =h'^{xo, Vo) + h]{xo, ?/o), for j = 1, 2, 3. 

Next, we treat one by one each of these terms. We shall assume everywhere, unless 
otherwise specified, that Xo,?/o > 0- To shorten the notation, we shall use the variables 
Zq and zq for either xq or yo, with the meaning explained in the Introduction for z,z and 
u,v. In the same way, we extend the notation for the characteristic functions (p,ip for 
conditions expressed in terms of zo,zo. 



6.1. The density h\{xQ,yo). In this case r = 1 and k = > 4 is even, although the 
first part of the calculation holds true more generally. Then, by (j221), we know that 

o 

Ck '■= C'k(l) = (^^O) ) , and by the definition of 7^ we know that Ck G if and only 
if the following conditions hold simultaneously: 



1 - ^ < 3;o < 1, 



3:0+1 ^ XQ+yo XQ+l 
k+1 ^ k ^ k ■ 



Since we assumed that xo,yo > 0, this translates into the equivalence 

xo + yo 



Ck & Tk 



< Xo, 2/0 < 1 and k < 



For the edges of 7^, we obtain 



Ckil)edTk\ViTk) 



and for the vertices of 7^., we have: 



1 - min(xo,?/o) 



for k>2. 



« - 5 < Xo< ^_^2' 

or 



k = ^f^, for j^<xo<l; 

l—yo ' k+2 " ' 

or 

xo = 1, for 1^ < 2/0 < 1; 



or 



yo = 1, 



Ck e V{Tk) 



'k — 1 
{xo,yo) G <! {-^^^ 



for 



k~l 
k+1 



< Xo < 1, 



k k 



A; + 2' fc + 2/' \ ' k + 1 



1,^ — -]\ (1,1) !> =Uk 



22 C. COBELI AND A. ZAHARESCU 

Then, using Table |21 for h^, the corresponding sum from the right-hand side of ()26j] gives 



hi{xo,yo) = ^ jvi^k < < xo,yo < l) 

k>4 
k even 

, 1 / /, xp + yo k-1 k \ 

fc>4 U I 

A; even 

+ ^(.0 = 1 if|^<^o<l)} (27) 

r 2k + i f k-i _ \ 

k even 

k + 2 ^( k \ 2^( \\ 



Here, for a given (xq, ^0)5 in each sum the number of nonzero terms is finite, at most equal 
to one in the last two sums that correspond to points on the border of W^. 



6.2. The density h\{xp^yp). Here r = 2, k = (1,/) or k = (A;,!), with A;,/ > 5, both 
odd. We assume first that k = (1,/). The center of Vi^i(2) has coordinates Ci^i : = 
Cu{2) = (xo, ^2^), by (j22D. Then Ci,/ G T^^i if and only if the following conditions hold 
simultaneously: 



2 ^ i-l ^ 



This gives 



G r^^i ^ { for / > 3. 

< / < min f 2o±M±l ^\ 

l-zo \^ 1-xo ' 1-yo / ' 
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Next, the conditions for the open edges of Ti i are 

(/ + l)xo + 2yo = l-l, for < Xq < 
or 



xo + lyo = I - 1, 
or 

Ixo + yo = I - 1, 
or 

yo = 1, 



for |^<xo<^; 



i-i. 



for < < ^; 



for L^<xo<^. 



For the vertices of Ti i, we have 



/ + 1 



/ + 2' / + 27' 
/ - 1 / - 1\ // - 2 



J + l'/ + l7' V / ' 

The symmetry allows us to use the zqj^o notation to collect the contribution of terms cor- 
responding to (1, /) and (/c, 1) in the same formula. Then, using Tabled the corresponding 
sum from the right-hand side of (^^1) gives 

2 / _ ^ 1+zq , . f2zQ + ZQ + l ^0 + 1' 

^ < / < mm I , 

Zq V 1 — 1 — ^0 



h2{xo,yo) = > 7 T^[zo,zo < 1; — < I < mm — , - 

-^/-IV I- zo \ 1- zo 1 

k odd 

+ Yl (7^1) {'^((^ + + 2zo = l-l, if ^ <zo< 



1-2 



l>5 
I odd 



1 + 2 



+ E 

/ odd 



2(/-l)(/ + l) 



, _ ^ - 2 _ 

+ (p[zo + lzo = l-l, if y-p^ < 2o < 1 



/-3 _ ^ 



(2^ 



+ 



2/2 + 5/ + 1 
4(/ - !)/(/ - 1)~ 



1-2 
1 + 2 



] Zq 



1 + 2 



l^f 21 + 1 ( 1-2 

+ tV'Uo = r— r + 7777 — ttV^^o = ; ^0 = 1 



I 



1 + 1) 4/(/-l) 
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Here, for a given (xq, ^0)5 in each sum the number of nonzero terms is finite, at most equal 
to one in the last two sums that correspond to points on the border of Ui^i^i. 



6.3. The density h^{xo, yo)- Now r = 3 and k = (1, /, 1), with / > 6 even. The center of 
n,/,i(3)isCi,z,i :=Ci,;,i(3) = (a:o, ^'~',y"° ),cf. ((221), and Tu,! = Ti,,. Then C^i G fi,,,i 
if and only if the following conditions hold simultaneously: 



(;_l) ('-i)^o+TO _i 



<^o< TTT 



XQ+l ^ (f-l)a:o+yo 
2 ^ 1-2 ^ 



This can be rewritten as 



a;o-2/o + 2 < /(I - yo), 

/(1-yo) <2(l + xo), 

/(l-a;o) < 2(2/0 + 1), 

l/o-a;o + 2 </(l-xo). 



The conditions for the open edges of T\^i^\ are 



/xo + 2yo = I - 2, 
or 

2xo + lyo = I - 2, 
or 



for |^<xo<j^ 



for fjj < a;o < — • 



(/ - l)a;o + yo = l-2, for ^ < xo < 



i-i. 



or 



^xo + (/ - l)?/o = 1-2, for < xo < 



For the vertices of 7i / 1, we have: 



(a;o,i/o) e 



/-3/-l\ /I -2 1-2 



l + l'l + lJ' V/ + 2'Z + 2/' 

/-I /-3\ d-2 1-2 
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Then, using Table |2l the corresponding sum from the right-hand side of fj26p gives 
h3{xo,yo) = j^^(^zo-zo + 2 < 1(1 -zq); 1{1 - zq) < 2(^o + 1)^ 
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Z>6 
k even 



{1-2) ^ ^{lzo + 2zo = 1-2, if <zo< 



l>6 
I even 



1-2 
1 + 2 



1-3 

TTT 



( w , .^-3 _ 1-2 

+ V9 ( zo + - l)-2o = ^ - 2, if r < -20 < 



/ + 1 



/ even 



2/ - 1 
4/(/ - 1 



/-3 



/ + 1 



^0 



/ - 1 

TTT 



+ 



/ + 2 
2Z(/ -2) 



/-2 



+ 



/ + 2/ 2(/-2)(/-l) 



1-2 



Uo 



(29) 



Here, for a given (xq, yo), in each sum the number of nonzero terms is finite, at most equal 
to one in the last two sums that correspond to points on the border of 

6.4. The contribution to the density of points of type T(r), r > 5. In this section 
we find the tail density, which we define to be 



r>5 



For any r > 5, there exist only two admissible r— tuples k, symmetric to one another. 
Let k = (1, 2, . . . , 2, 3) be such an r— tuple. Then, since p.r(k) = 2 and pr-i{2, . . . , 3) = 
2r - 1, relation ^ produces := C^{3) = {xq, i2Lzl^2±my Then Ck G Tk if and only 
if the following conditions hold simultaneously: 

2 < Xo < 



2r-l 



2r-3 



This can be rewritten as 



yo < 1, 

1 < 2xo + yo, 

2 < (2r + l)xo + yo, 
2>{2r- l)xo + yo- 



26 C. COBELI AND A. ZAHARESCU 

The conditions for the open edges of 7^ are 



Ck e 9Tk \ V(Tk) 



(2r + l)xo + 1/0 = 2, for ^ < Xq < 



2r-l . 

i+2 ' 



or 

+ Z/o = 1, 
or 



for 273^ < a;o < 2733; 



(2r - l)xo + j/o = 2, for ^^i^ < a;o < 2^1: 



or 



for tt::^ < a;o < ^ 



2r+l 



2r-l" 



For the vertices of 7k, we have 



1 2r - 3 



2r + l' 7' V2r-l'2r-l/' 
1 2r-5^ 



Wk. 



.2r - 3' 2r - 3/' V2r - 1' 
We use the notation with the variables zo,Zq and the symmetry to write in the same 
formula the contribution to 5'(xo,?/o) of all the terms corresponding to the r— tuples 
(1, 2, ... , 2, 3) and (3, 2, ... , 2, 1), for r > 5. Using the information from Table El we 
obtain 

g'^ixo, yo) = ^'f{zo<l< 2zo + Zq; (2r - l)zo + ^0 < 2 < (2r + 1)2:0 + ^0) 



r>5 



<Zo< 



2r - 1 



+ J]iL((2r + l>o + ^o = 2, if^ 

+ ^(2zo + zo = l, if <:z^<:-l 
+ ^((2r - 1)^0 + Zo = 2, if < zo < 



(30) 



8(2r + l)^V° 2r + 
14r + 9 



16(2r + 1) 



1 



2r - 3 



2r -3 

+ 7T^[Zo 



1 



Zo 



2r - 1 2r - 1 

2r - 5 



16(2r - l)'^ V " 2r - 3 2r - 3 



+ 



4r - 1 
8(2r - 1) 



2r - 1' 



^0 = 1 
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Here, in each sum, for a given {xo,yo), at most one term is nonzero. 
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6.5. The baby puzzle. Let i > 5 be a positive integer. Then, we observe that the 
polygons Ui^i, Ui^i, fit perfectly into the quadrilateral with vertices 

2 — 3 \ — 1 i — 1\ / i — 3 



(see Figure EJ. Another nice aspect of this matching is due to the fact that — 1) = 
P2(l, i) = P2{'i, 1) = ^3(1, ^ + 1, 1), so 0i may be viewed as region with constant density at 
level i, say. Let us notice that D (257 D (Sg D, . . . , so gi'^'', "w) is not locally constant on 
(5i, but the support of g{u, v) is a superposition of quadrilateral steps of constant density 
(provided we show that a similar property holds for the remaining polygons Wk). Putting 
together d^Zj), (j2Hl), (EHl), we obtain the next result. 

Proposition 1. For any (xo,yo) ^ [0? have 
h'^ixo, yo) ■=hl{xQ, yo) + h^ixo, yo) + hl{xo, yo) = 

i>5 ^ 
i odd 

+ 4„=1. if <.„<!)} (31) 



j>5 
i odd 



+ E| 2(.-i)(. + i) K"° = 7TT'^° = ' 



i>5 
j odd 



« + 3 « - 1\ 1 ~, 

+ 2(z - l)(z + = TT^J + 2(z3T)'^(^o = 1 



Proof. We only need to check the equality at the matching corners. We have: 

^"(7^' =2("^'*'^^'^'' +"i.i+i.i(fTi' 1+2)) 

_ 1 / + 5i + 1 2i + 1 \ _ 1 
~2 V2(i-l)i(i + l) ^ 2i{i + 1) y ~ i - 1 ' 
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and 

_ 1 / 2i + 1 2z - 1 \ _ 1 
~2 \2i{i - 1) ^ 2i(i - 1) ) ~ i - 1 ' 

equal each to half of the interior density, since they are on the open edges of C5j. By 
symmetry, we have the same result at + 2), (z — 2)/(i + 2)) and at (l, (z — 2)/-i). At 
the interior matching point, we have 

~2 Vi(i - 1) ^ 7 ^ i(i - 1)7 ~ i - 1 ' 
which concludes the proof of the proposition. □ 

6.6. The big puzzle. We group the terms of lower orders into 

/(xo, yo) = h'^ixo, yo) + h'^i^xo, yo) + g^^xo, yo), 

where h'^{xo,yo) = hf{xo,yo) + h'l{xo,yo) + h'l{xo,yo). It remains to find h'^{xo,yo) + 
g'i{xo,yo), as /i"(xo,?/o) was the object of Section In fact, the calculations are special 
cases of those already performed in sections 16. 1116. 41 Here, we shall show that the sum of 
h'^, and can be combined into a simpler formula, similar to that of h^. 

hetM = {{2)- (1,3); (3,1); (1,2,3); (3,2,1); (1,4,1); (1,2,4,1); (1,4,2,1); (1,2,2,3); 
(3,2,2,1); (1,2,2,2,3); (3,2,2,2,1);...}. The key point in the matching that occurs 
among the supports of different components of /i*^ + (?4 + g^ is the fact that the kernel 
is constant, equal to 1, for any k contributing to the sum. This follows by the equality 
Pr(k) = 2, for any k G M. 

Let be the quadrilateral with vertices (0,1); (1/3,1/3); (1,0); (1,1). It turns out 
that is the support of h"^ + g4^ + g^, and it looks like a mosaic (that is, no superpositions 
over interior points occur) composed by all polygons Wk, with k G M.. Perfect matchings 
occur (see Figure IH)) getting particular cases of the general relations, as follows. (In 
the nonsymmetric cases we give the statement only for the polygons situated above the 
first diagonal.) The polygons ?7i,2,2,3 and t/1,2,3 U Wi_2,4,i are given also by the formula 
for Wi_2,...,2,3! with r = 4 and r = 3, respectively. More significantly, on the one hand 
^1,4,1 = ^1,2,4,1 U Ui^i^i U Wi_4_2,i is given by the formula for Wi^^^i, with I = 4 and, on the 
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other hand, Ui^z = ^1,3 U Wi,2,3 U Wi_2,2,3 U ^1,2,2,2,3 ■ ■ ■ is given by the formula for Ui^i, with 
/ = 3. Then © is composed by Wi^s, Wi,4,i, ^3,1, and Ui in the same way as 
Wj^i, and completed the baby puzzle, 0j. Next we summarize the contribution to 
g{xQ,yo) of all tuples k G A^. 

Proposition 2. For an?/ (xo,?/o) ^ [0, 1]^, 'we /iave 
/"(a;o,Z/o) : = /i'^lxcZ/o) + fl'4(a;o, Z/o) + g''{.xo,yo) = 

= '^[zo<l; 1< 22:0 + ^0) + 

1 1 (32) 

+ -<^(22;o + zo = l, if < zo < 1/3) + -<^(-2o = 1, if < < l) 

3 3 _ 1 _ 

+ T^^(^o = 0; Zo = l) + -(p{zo = 1/3) + -(p{zo = 1) . 
lo 8 4 

Proof. We just check the equality at the matching vertices using entries from Table El 

Beginning with the border of 0, for points on the top edge, we have: 

/"(1/3, 1) =^(«2(l/3, 2/3) + ai3(l/3, 1)) = ^ + ^) = ^ ! 
/"(1/5, 1) =^(«i3(l/5, 4/5) + ai23(l/5, 1)) = ^ + ^) = ^ 5 
9'^{l/7, 1) =^(«i23(l/7, 6/7) + ai223(l/7, 1)) = ^ (i + i) = ^ ^ 

9'^''i2^^^) =2 ("^^■■■^^^2;^' 2^'' +"122-23(2^,1)) 

1 / 4r + 1 4(r + 1) - 1 \ 1 ^ 
- ' + — ^- = -, for r > 4. 



2V4(2r + l) 4(2(r + l)-l)y 2 
On the edge with endpoints (0, 1), (1/3, 1/3), we have: 

g "(1/5, 3/5) =- («i23(l/5, 4/5) + «i223 (1/5,1) (1/5,4/5)) 





'13 


5 




1 




.21 ^ 


^6^ 


^24) 


2 



/"(1/7, 5/7) =i («i223(l/7, 6/7) + ai2223(l/7, 1)) = ^ + ^) = ^ ^ 

^'^"(i;^!, izf) =2 ("12-23 (2^31,^^) +ai22-23 (2731,1)^ 

1 / 14r + 9 2(r + 1) - 3 \ 1 

+ r = -, for r > 5. 



2V8(2r + l) 8(2(r + l)-l)7 2' 
These are all equal to half of the kernel, as expected. The same results hold for points 
symmetric with respect to the first diagonal. 
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At the vertices of 0, we have: 

=-a2(l,l) = -■- = -] 

/"(1/3, 1/3) =^ (ai4i(l/3, 2/3) + ai24i(l/3, 1) + ai42i(l/3, 2/3)) 

_1 /3 3 3 \ _ 3 
~2V5^40^40y "18' 

At the corners (0, 1) and (1, 0) the problem is a httle bit more comphcated. Here, as the 

inside hmit in (fT^ is taken over Q, the parallelogram V]<i{r) will cover completely 7k, for 

infinitely many k from the sequences with 2's embraced by 1 and 3, while some of these 

Tk's will be covered partially. Though, we have a quick shot solution to the problem of 

finding ^''^"(0, 1) and 5''^"(1,0) using, on a larger scale, the property of the parallelogram 

used in for the quadrilateral with vertices: (0, 1); (1/3, 1/3); (1/2, 1/2); (1/3, 1) and 

the fact that in C5 the index has everywhere the same value 2. We obtain 

^7(0, 1) =^ [2 - iai4i(l/2, 1) - ai3(l/2, 1) - ai3(l/3, 1) 

-ai24i(l/3,l)-^ai4i(l/3,2/3) 

_1 / _ 1 2 _ 1 _ 7 3 1 3\ _ 3 
~2\ ~ 2 ' 3 ~ 3 ^ 12 ~ 40 ~ 2 ' Sy 16 ■ 

Finally, we complete the puzzle with Ti_4 1, the last piece. We have 

/"(1/2, 1/2) =^(a2(l/2, 1/2) + ai3(l/2, 1) + a3i(l/2, 1/2) + ai4i(l/2, 1)) 

_l/2 1 1 2 
~2 V3 ^ 3 ^ 3 ^ 3 

/"(2/7, 4/7) =^(ai3(2/7, 5/7) + ai23(2/7, 1) + ai4i(2/7, 5/7) + ai24i(2/7, 1)) 

_1 /19 11 11 19\ _ 
~2 V30 ^30^30^307 ' 

and, by symmetry, /"(4/7,2/7) = 1. We see that g'^'' takes at (1/2, 1/2), (2/7,4/7) and 
(4/7, 2/7) the value 1, equal to the kernel, as needed, since these are interior points. This 
completes the proof of the proposition. 

□ 

6.7. Completion of the proof of Theorem [H We first notice that can be obtained 
as a particular case of the formula for C5j, with i = 3. Moreover, the right-hand side 
of can also be obtained if we put i = 3 into the generic terms of the sums on 
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the right-hand side of (jHTj) . Thus, we employ Propositions ^ and |21 to obtain the sum 
gi^o, yo) = g'^^i^o, yo) + h^i^o, yo), which completes the proof of Theorem [TJ 



7. Proof of Theorem [21 

Let I C [0, 1] be fixed and let a'/q',a"/q" be consecutive fractions in Let us first 
see how one translates the condition a'/q' G X in terms of the two variables q', q". Since 
a"q' — a'q" = 1, it follows that a' = —q" (mod q'). Here q" is the representative from 
[0, g' — 1] of the inverse of q" modulo q' . Then, one immediately derives that a'/q' G X 
if and only if q' — q" G q'T. Similarly, we get that a" /q" G X if and only if q' G g"X, but 
here the inverse is taken modulo q" . We remark that one of these conditions is almost 
redundant, since a'/g' G X ensures that a" /q" G X, also, except for at most one pair g', g", 
and conversely. Then, imposing only one of these two conditions, one may neglect this at 
most one term in the corresponding computations below, and absorb it in the error term. 

The proof of Theorem |21 follows the same steps from the beginning of the proof of 
Theorem d We have to find the ratio of the number of elements in the set 



. 1 < g', q" < g, gcd(g', q") = 1, g' + q" > Q, q' mod q" G g"X 
Bl{r)={{q',q')eN': 



q' even, q" odd; k(g',g") G A{r), {q',q^{r)) G Q ■ □ 



and the cardinality of i^^^^^^^^- Then the turning point is the analogue of relation 
which becomes 



Qj. (x, y) dxdy = lim . (33) 



Q,cvcn 



Q ,cvcn 



This will complete the proof, provided we show that ^^Bqlr) ~ |X| ■ ^Bqi^r) and 
|X| ■ as Q ^ oo. One should observe that, with the notations from Sectional we 

haveSQ(r) =i3g''^(r). 

To proceed, we estimate i?^ 

Lemma 4. For any subinterval 2 C [0, 1], we have 
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Proof. The cardinality of can be written as 



q even, q odd, q' mod q G g X 
G (<5 — Q', Q] : gcd(x, g) = 1, x even, x mod q G gX} . 

l<q<Q 
q odd 



(35) 



Notice that the condition q' G g"X introduces randomness in the positioning of points 
from the right hand side of (|35|) . In order to estimate the terms added in the sum, we 
write them using exponential sums, separate the main term, and employ classical bounds 
for Kloosterman sums (see ^H], Thus, they can be written as 

Q-q<x<Q y&qX ^ k=i ^ 
gcd(a:,q)=l 
X even 



^ ^ k=l yeql ^ Q-q<x<Q ^ 

gcd{x,q) = l 
X even 

= ^. f., HI + 

q 2 ^ ^ 

in which 0(g) is Euler's totient function. Then, substituting in (j35j) . we obtain 



l<g<Q 
q odd 



It remains to estimate the sum from (j36|) . This is 



E '-f-^- E E^-^-E^E^-^. 

l<g<Q ^ l<q<Q d\q d=l qi=l 

q odd q odd d odd odd 

" M<i) 1 , (37) 



i odd 

2Q' 



^2 



+ 0(QlogQ), 
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since, via the Euler product, we find that 



d=l p>3 
d odd 



Now, the lemma follows by inserting the estimation ()37|1 on the right-hand side of □ 

We remark that the size of the error term in (jH^ may be slightly lowered, but this is 
not essential for our needs. 

Next, for any ^7 C and X C [0, 1], we denote 

^e^en,odd(^) := # { (a^, Z/) ^ 1^ fl : x even, y odd, gcd(a;, y) = 1, x mod y G yl] . 
Then, in the spirit of Lemma EJ we get, more generally, the following result. 
Lemma 5. Let R > and Q C [0, R] x [0, R] be a convex domain. Then, we have 

NLn,oM = 1^1 ■ ^even,odd(^^) + 0(i?3/2+^) _ 

Proof. The proof follows the first part of the proof of Lemma IH We have: 



Arx /ON // ) / N o^^t2 gcd(gi,g2) = 1, gi even, odd, 

Cen,odd(^) 92) G^^nN^: _ 

I qi mod q2 G q2l 

= ^ G J(g) : gcd(x, q) = 1, x even, x mod q G gX} , 

l<q<R 
q odd 



(38) 



in which 3{q) := QCl {y = q}. Evaluating the size of the terms here, using again the same 
estimates for Kloosterman sums, we find that they are equal to 

gcd(a;,g)=l 
X even 

The required result follows by (jHHj) and (jSHI), and the fact that 



iVe.en,odd(n) = ■ + 0(i?) . 

i<g<R 

q odd 



□ 
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Now we have all the tools needed to complete the proof of Theorem |21 Since ^^^^ = 
QV^' + O(QlogQ), by Lemmallwe find that = |J| ■ + 0(Q3/2+.). Qn 

the other hand, using notations from Sectional we find that B^{r) = -^^en odd(^Q('^))' 
where ^gir) = QQ{xo,yo,v){r) is given by 

{„ l<x,y<Q, x + y>Q, 
k(x, y) G A{r), {x, x^{x, y)) e Q ■ D 

The set ^qir) is in general not convex, but it is a finite union of boundedly many convex 
sets, as Q ^ oo, the number of these convex sets depending on the given point (a;o, yo) ^ 
[0, 1] X [0, 1]. Then, by Lemma it follows that 



#i3g(r) = |T|-#i3Q(r) + 0(.„,,)(g='/2+ 
concluding the proof of Theorem |21 
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